Systematic study of 91011 Li with the tensor and pairing correlations 
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We make a systematic study of Li isotopes (A = 9, 10, 11) in the tensor optimized shell model 
for 9 Li and treat the additional valence neutrons in the cluster model approach by taking into 
account the Pauli-blocking effect caused by the tensor and pairing correlations. We describe the 
tensor correlations in 9 Li fully in the tensor-optimized shell model, where the variation of the size 
parameters of the single particle orbits is essential for getting strong tensor correlations. We have 
shown in our previous study that in 10,11 Li the tensor and pairing correlations in 9 Li are Pauli- 
blocked by additional valence neutrons, which make the p-shell configurations pushed up in energy. 
As a result, the s 2 valence neutron component increases to reveal the halo structure of 11 Li and the 
inversion phenomenon of the single particle spectrum in 10 Li arises. Following the previous study, 
we demonstrate the reliability of our framework by performing a detailed systematic analysis of the 
structures of 9,10,11 Li, such as the charge radius, the spatial correlation of halo neutrons of n Li 
and the electromagnetic properties of Li isotopes. The detailed effects of the Pauli-blocking on the 
spectroscopic properties of 10 ' n Li are also discussed. It is found that the blocking acts strongly 
for the 11 Li ground state rather than for 10 Li and for the dipole excited states of n Li, which is 
mainly caused by the interplay between the tensor correlation in 9 Li and the halo neutrons. The 
results obtained in these analyses clearly show that the inert core assumption of 9 Li is not realistic 
to explain the anomalous structures observed in 10 ' n Li. For the dipole excitation spectrum of n Li, 
the effect of the final state interactions is discussed in terms of the dipole strength function. 

PACS numbers: 21.60.Gx, 21.10.Pc, 27.20. +n 



I. INTRODUCTION 

The halo structure in 11 Li was completely outside of 
the common sense of nuclear physics |l[ . Because of this 
reason, many experiments have been performed on n Li 
and the surrounding nuclei in order to understand the 
physics causing this interesting phenomenon. We sum- 
marize here the available experimental data accumulated 
for 11 Li and surrounding nuclei until now: a) The halo 
neutrons have almost equal amount of the s-wave compo- 
nent with respect to the p-wave component Q . b) More 
precise measurements and analyses of r.m.s. radius are 
obtained |3|, |4|. c) The El strength distribution has a 
large peak near the threshold @, HJ0, 1 • d) The charge 
radius is larger than that of 9 Li [9|, e) The elec- 

tromagnetic pro pert ies are observed to reveal the exotic 
structure of 11 Li[lj. 

The challenge for the theoretical work is the large s- 
wave component for the halo neutrons, which means that 
the shell gap at N = 8 has to disappear by some reason. 
However, the mean field treatment of a central force is 
not able to provide the disappearance of the N = 8 shell 
gap, because the N — 8 magic number is produced by 
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the shell structure due to the central potential. Hence, 
there have been many theoretical studies performed for 

but con- 
ventionally, most of them had to assume that the lsi/2 
single particle state is brought down by hand to the Opi/2 
state in order to match with the experimental situation. 
It is therefore the real challenge for theorists to under- 
stand this disappearance of the N = 8 shell gap, to be 
called the s-p shell gap problem. 

The halo structure of 11 Li should be related with the 
states involving the Is state (Is states) and those involv- 
ing the Op-state (Op states) in 10 Li. Several experiments 
suggest that non-resonant, virtual ls-states are located 
close to the 9 Li+n threshold energy together with the 
nearby but higher Op states [24|, |25|, l26[ . The formation 
of the low-lying virtual ls-state in 10 Li is definitely im- 
portant to understand the inversion phenomenon of the 
Is and Op states in the N — 7 isotones. However, this 
inversion phenomenon is still an unsolved problem po| 
in addition to the halo structures of 11 Li. Hence, we 
need an important mechanism to explain the inversion 
phenomenon of the s-p states and the virtual ls-state in 
10 Li consistently with the structures of n Li. 

It is well known by now that the tensor interaction 
plays an important role for light nuclei. The deuteron is 
bound due to the strong coupling between the s-state and 
the d-state due to the tensor interaction. In 4 He, the ten- 
sor interaction produces a strong coupling of the s-state 
and the d-state in the relative wave function and even- 
tually the contribution of the tensor interaction becomes 
energetically comparable to or even more than that of the 
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central interaction [27], [28| . Very recently, the Argonne- 
Illinois group made a thorough study for light nuclei until 
mass number of about A=10 using the Green's function 
Monte Carlo method in the few-body framework. They 
take the bare NN interaction including the tensor inter- 
action with a small but important three body interaction 
for quantitative reproduction of light nuclei [29| . They 
found not only very good reproduction of binding ener- 
gies, but also an important fact that the pion exchange 
interaction (predominantly tensor interaction) provides 
about 70^80% of the two-body attraction for those light 
nuclei. Hence, a proper treatment of the tensor interac- 
tion seems essential for the understanding of the nuclear 
structure. 

In our recent study [13, HH , we have developed a the- 
oretical framework of the tensor-optimized shell model 
(TOSM) to treat the tensor interaction in the shell model 
basis explicitly. This finding is very important since we 
now have a tool to systematically understand the role 
of the tensor interaction and its correlation in the shell 
model language. The first place to use this method was 
the 4 He+n system, where we are able to show that the 
tensor correlation provides a considerable amount of the 
P1/2-P3/2 splitting energy in 5 He [13, H2, H3, HIl • Hence, 
we can apply this method for the present shell gap prob- 
lem. In our previous paper [1||, we presented our basic 
idea and the formulation to treat the tensor correlation 
in the 9 Li+n+n system. In the present paper, we shall 
demonstrate the reliability and the validity of this frame- 
work for Li isotopes by showing various physical quanti- 
ties observed not only for 11 Li, but also for the neighbor- 
ing isotopes 9 Li and 10 Li. We have applied TOSM to the 
9 Li core to describe the tensor correlation explicitly. It 
is found that the (0si/2) _2 (0pi/2) 2 of the 2p-2h excita- 
tion of a proton-neutron pair has a special importance in 
describing the relative sd coupling of the tensor interac- 
tion and the tensor correlation in 9 Li, which is the same 
result as obtained for 4 He [H H, Ha, HI] . For n Li, we 
solved the coupled 9 Li+n+n three-body problem with 
the tensor correlation in 9 Li 35]. The tensor correlation 
arising from the 2p-2h excitations in 9 Li is Pauli-blocked 
by additional two neutrons, which makes the p-shell con- 
figurations pushed up in energy. As a result, the s-p shell 
gap becomes prominently small and the s 2 component for 
the last neutrons increases largely around 25%. 

The tensor correlation is, however, not enough for 
quantitative description of the s-p shell gap problem. We 
have to add further the paring correlation among the p- 
shell states. Hence, we have solved a coupled problem 
with the tensor and pairing correlations in 9 Li. We con- 
firmed then that the neutron pairing correlation in 9 Li 
works coherently with the tensor correlation and pro- 
duces the similar blocking effect in 11 Li as the case of the 
tensor correlation [13, H3, 53] ■ As the final result with 
the tensor and pairing correlations, the s-p shell gap dis- 
appears and the s 2 component for last neutrons naturally 
increases as around 50% which is enough to produce the 
halo structure of 11 Li. Our theory further reproduces the 



large charge radius and the Coulomb breakup strength of 
n Li, and the inversion phenomena of s-p states in 1 Li. 

So far, there have been appearing a number of exper- 
imental data on 11 Li and neighboring nuclei. Hence, in 
this paper, we shall further perform a detailed and sys- 
tematic analysis of the structures of the Li isotopes on 
the basis of the tensor optimized shell model (35[. We 
investigate the Pauli-blocking effect on various physical 
quantities of 10 ' n Li, such as the 10 Li spectra of s and p 
states and the ground state properties of 11 Li including 
the spatial correlation of halo neutrons. The charge ra- 
dius, the electromagnetic properties of 9 Li and n Li are 
also investigated. For the three-body Coulomb breakup 
of n Li, we confirm the low-energy peak in the strength 
and the origin of this peak is investigated from the view 
point of the final state interactions. The difference of the 
Pauli-blocking effects in 10 Li and n Li is also discussed, 
which is caused by the coupling of the 9 Li core excitation 
and the halo neutrons. 

This paper is arranged as follows. In Sec. [ill wc explain 
the wave functions of 9 ' 10 ' n Li, which describe the tensor 
and pairing correlations explicitly. In Sec. IIII1 based on 
the obtained results of 9 Li, we show the role of the Pauli- 
blocking effects coming from the tensor and pairing cor- 
relations on 10 ' n Li and show various properties of these 
nuclei. A Summary is given in Sec. HVl 

II. MODEL 

A. The tensor optimized shell model for 9 Li 

We first explain the tensor optimized shell model 
(TOSM) for 9 Li The Hamiltonian for 9 Li is given 

as, 

9 

ff( 9 Li) - » '''••»<•' ' W 

i=l i<j 

Here, ti, to, and are the kinetic energy of each nu- 
cleon, the center-of-mass (cm.) term and the two-body 
NN interaction consisting of central, spin-orbit, ten- 
sor and Coulomb terms, respectively. The wave func- 
tion of the 3/2~ state in 9 Li is described in the tensor- 
optimized shell model [13, [HI- We express 9 Li by a multi- 
configuration, 

N 

*( 9 Li)=5><I>f ~, (2) 

i 

where we take the Op-Oh state of 9 Li as (0s) 4 (0p 3 / 2 ) 5 
and consider up to the 2p-2h excitations from the Li 

3 /2 — 

core state configuration within the Op shell for $/ in 
a shell model type wave function, and N is the configu- 
ration number. Based on the previous study [13, HH Ha] , 
we adopt the spatially modified Gaussian function as a 
single-particle orbit and treat the length parameter b a of 
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each orbit a (Os, Opi/2 and 0^3/2) as variational param- 
eters. 

In the tensor-optimized shell model, we solve the fol- 
lowing variational equations with respect to the total 
wave function ^( 9 Li) as, 

d{*iff( 9 Li)-£|*) = o d{m( 9 ^)-m) = (3 

dcii db a 

Here, E is the total energy of 9 Li. We determine {oj} 
and the length parameters {b a } of three orbits. We 
can optimize the radial form of single-particle orbits ap- 
propriately in the wave function in Eq. ([5]), such as 
to describe the spatial shrinkage of the particle state, 
which is important in realizing the tensor correlation 
[H M, HE m, H [H, EO]. We add further the neutron 
paring correlation for the p-shell orbits. 

In the description of the tensor correlation, in princi- 
ple, we can work out a large space to include the full 
space effect of the tensor force by taking 2p-2h states 
with very high angular momenta [3l| . In order to avoid 
large computational efforts without loss of the physical 
importance in the result, we restrict the 2p-2h shell model 
states within the p-wave states for the description of 9 Li 
with the single Gaussian basis. Instead, we increase the 
tensor matrix elements by 50%. The justification of this 
enhancement factor is explained as follows, which has 
also been discussed in Ref.[31]. In TOSM, addition of 
the higher configurations beyond the p-shell increases 
the total amount of the 2p-2h components in the wave 
function, while these configurations share the total 2p-2h 
components and reduce the 2p-2h components includ- 
ing the p-shell configurations [3ll| . On the other hand, 
the spatial behaviour of each single particle basis can 
be improved also using the Gaussian expansion method 
(GEM), from the compact Gaussian form to more suit- 
able one. This GEM effect was found to enhance each 
2p-2h component [3ll. [35| . Both effects can be included in 
the wave function of the core in TOSM. When the Pauli- 
blocking effect is considered by adding an extra neutron 
into the particular occupied orbit of neutron in the core, 
the mixing of higher configurations reduces this block- 
ing because of the sharing of the 2p-2h components. On 
the other hand, the GEM effect enhances the blocking, 
because the Pauli-blocking is proportional to the spatial 
overlap of the neutron wave functions between the core 
and the extra neutron. Actually, it was checked that in 
the scattering problem of the 4 He+n system, the com- 
bined effect of the higher configurations and GEM on 
the Pauli-blocking is simulated using the enhanced ten- 
sor matrix elements with a single Gaussian basis [32l |35|. 
Therefore, in the present study, we adopt this approxi- 
mation of the tensor correlation for 9 Li. 

We discuss here the interaction for 9 Li; in iJ( 9 Li), 
where we use the shell model wave functions for the 9 Li 
core in Eq. ([2]). Since our main interest in this work is 
to investigate the role of the tensor interaction on the 
two-neutron halo formation, we describe the tensor cor- 
relation in addition to the pairing correlation in 9 Li. We 



use the GA interaction proposed by Akaishi [3(| 5H 
for Vij in Eqs. ([1]), ([8]) and ((4]). This effective interac- 
tion GA has the tensor force obtained from the G-matrix 
calculation using the AV8' bare force keeping a large mo- 
mentum space |42| . [43| . Actually, it was shown in the 
previous paper [3l| that the tensor force matrix element 
using GA for the relative sd coupling almost reproduces 
the results of AV8'. This means that we employ the ef- 
fective tensor force in GA, which almost has the same 
strength of the bare tensor force. In GA, the obtained 9 Li 
wave function in Eq. ^ shows smaller matter radius than 
the observed one due to the high momentum component 
produced by the tensor correlation [3(| 0, Hence, 
we phenomenologically adjust the central force, which is 
done by changing the second range of the central force by 
reducing the strength by 21.5% and increasing the range 
by 0.185 fm to reproduce the observed bindin g en ergy 
(45.3 MeV) and the matter radius of 9 Li[3(J [Hill]. 

B. 10 Li as 9 Li+n system in TOSM 

For 10 Li, the Hamiltonian for 9 Li+n is given as 

H( 10 U) = H{ 9 U) + T + Ti - Tg> + V cn , (4) 

where H( 9 Li), T t and T^ 2) are the internal Hamiltonian 
of 9 Li given by Eq. (fl]) , the kinetic energies of each cluster 
(T for 9 Li, the same as to in Eq. |T])), the cm. term of 
two cluster systems, respectively. V cn is the 9 Li core-n 
interaction. The wave function of 10 Li with the spin J is 
given as 

N 

^( 10 Li) = Y. A {^' 2 \xi°{n)} J ). (5) 

i 

We use the 9 Li wave functions in TOSM and obtain cou- 
pled differential equations for the neutron wave function, 
X l7 °(n) in order to describe the relative motion between 
9 Li and neutron, where Jo is the total angular momen- 
tum of the last neutron part of 10 Li. To obtain the total 
wave function, ^ J ( 10 Li), we use the orthogonality con- 
dition model (OCM) [lj, HO, H] to treat the antisym- 
metrization between the last neutron and neutrons in 9 Li. 
In OCM, the neutron wave functions x are imposed to 
be orthogonal to the occupied orbits by neutrons in 9 Li, 

which depends on the configuration ( E ) ^ 2 in Eq. @. We 
obtain the following coupled Schrodinger equations with 
OCM for the set of the wave function, {Xi°( n )}i where 
' 1.----.V: 

N 

Yl [( T rel + V ™ + ^) % + M^) 

i=i 

x X j°(n) = E X i°{n), (6) 
A 2 = A IV'aXM (7) 

ae*i( 9 Li) 
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where /iy( 9 Li) 



netic energy of the relative motion for 10 Li. Aj is the pro- 
jection operator to remove the Pauli forbidden states tjjg 
from the relative wave functions of last neutron (iol |4q|. 
where ip a is the occupied single particle wave function of 
the orbit a in 9 Li. This Aj depends on the configuration 

$^/ 2 of 9 Li and plays an essential role to produce the 
Pauli-blocking in 10 Li. The value of A is taken as large 
as 10 6 MeV in the present calculation in order to project 
out the components of the Pauli forbidden states into an 
unphysical energy region. Here, we keep the length pa- 
rameters b a of the single particle wave functions as those 
obtained for 9 Li. The detailed treatment of the orthogo- 
nality condition was explained in Ref. [35[. 

The 9 Li-n potential, V cn , in Eq. ^ is given by 
folding an effective interaction, the MHN G-matrix 
interactionpjl, [H Hf| |4?], H^]. Any state-dependence is 
not introduced in the 9 Li-n potential. We introduce one 
parameter, 5, which is the second-range strength of the 
MHN interaction to describe the starting energy depen- 
dence dominantly coming from the tensor interaction in 
the G-matrix calculation [45|, |H, |4i| . In the present cal- 
culation, we chose this S parameter 0.1745 to reproduce 
the two-neutron separation energy S2n of n Li as 0.31 
MeV, since experimental situation is not so clear for the 
10 Li spectra due to the unbound system. 



is the ki- 



C. n Li as 9 Li+n+n system in TOSM 

For 11 Li, the Hamiltonian of the 9 Li+n+n three- 
cluster system, is given as 



2 2 

if ("Li) = H ( 9 Li) + ]T T k - + ]T W 

k=Q k=l 



(8) 



(3) 

where Tq is the kinetic energy of the c. m. term of 
three cluster systems. V cn ,k are the 9 Li core-n interaction 
(k = 1, 2) and V nn is the interaction between the last two 
neutrons. The wave functions of 11 Li with the spin J', 
are given as 

N 

v^'C'Li) = E- A {[** /2 "'X?( n »)] J '}- ( 9 ) 

i 

We obtain coupled differential equations for the neutron 
wave functions x J ° ( nn ) > where J' a is the total angular 
momenta of the last neutrons. To obtain the total wave 
function ^ J ( n Li), we use the orthogonality condition 
model (OCM) pjl, Gil Ell to treat the antisymmetrization 
between last neutrons and neutrons in 9 Li. In OCM, the 
neutron wave functions x are imposed to be orthogonal 
to the occupied orbits by neutrons in 9 Li, which depends 

on the configuration in Eq. @. We obtain the 

following coupled Schrodinger equations with OCM for 



the set of the wave functions {Xi°( nn )} f° r n Li, where 



N 

E 

3=1 



Sij + hij ( Li) 



k=l 



x Xj°( nn ) — E Xi°(nn), 



(10) 



where /ijj( 9 Li) = ($ 



3/2' 



\H( 9 Li)\$ 3 / 2 ). T^i is the total 



kinetic energies consisting of the relative motions for n Li. 
A^fc is the projection operator for each last neutron with 
an index k. 

We describe the two neutron wave functions x( nn ) m 
Eq. (JTOJ) precisely adopting a few-body approach of the 
hybrid-TV model [H III : 



X J i°{nn) 



(11) 



where £ v and £ T are V-type and T-type coordinate sets 
of a three-body system, respectively. The radial part of 
the relative wave function is expanded with a finite num- 
ber of Gaussian basis functions centered at the origin, 
which can describe the loosely binding wave function of 
neutron halos[4o|. 

Here, we discuss the coupling between 9 Li and the last 
neutrons, whose details were already explained in the pre- 
vious paper [nl [35[. We consider the n Li case. Asymp- 
totically, when the last two neutrons are independently 
far away from 9 Li (£y T — > oo), any coupling between 
9 Li and two neutrons disappears and 9 Li becomes its iso- 
lated ground state. Namely, the mixing coefficients {a,} 
are the same as those obtained in Eq. @: 



$ J '( n Li) 



*( 9 Li), X J °(nn) 



J' 



f( 9 Li) = 5>$ 



3/2" 



(12) 
(13) 



where {Xi° ( nn )} m Eq. © satisfies the following asymp- 
totic forms for i = 1, • • • , N 



xi°{nn) 



di ■ x J °(nn). 



(14) 



This equation implies that the asymptotic wave function 

of two neutrons Xi"i nn ) is decomposed into the inter- 
nal amplitude a, of 9 Li and the relative wave function 
X J °{nn). Equations (|T2]) - (p~4|) give the boundary condi- 
tion of the present coupled three-body model of 11 Li. 
Contrastingly, when the two neutrons are close to 9 Li, 
two neutrons dynamically couple to the configuration 
^3/2 o £ g-y sa (.j g £y m g Pauli-principle. This cou- 
pling depends on the relative distance between 9 Li and 
the neutrons and changes {a^} of 9 Li from those of the 
9 Li ground state. The structure change of 9 Li inside n Li 
is determined variationally to minimize the energy of the 
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n Li ground state. The dynamical effect of the coupling 
arising from the Pauli-blocking is explained in the results 
in detail. 

For the potential V nn between the last two neutrons, 
we take a realistic interaction AV8' in Eq. (JSj> . Our in- 
terest is to see the n-n correlation in the two-neutron 
halo structure. Hence we solve two-neutron wave func- 
tion x( nn ) with no assumption in the hybrid- TV model 
in Eq. (fTTj) . which can describe the short range correla- 
tion under the realistic interaction AV8'. 



III. RESULTS 
A. 9 Li 

We show first the results of 9 Li in TOSM, which is 
the basis for heavier Li isotopes. As the first step, we 
search for the energy minima of 9 Li with respect to the 
length parameters {b a } by solving Eq. Fig. [T]shows 
the result of the energy surface with respect to the 6o Pl/2 , 
where &os arid 6op 3/2 are already optimized as 1.45 fm and 
1.8 fm, respectively. We found two energy minima, (a) 
having a small value (0.85 fm) of bo Pl/2 and (b), which 
is a rather shallow local minimum, having a large value 
of bo Pl/2 as 1.8 fm, the same length of bo P3/2 . The con- 
tribution of the tensor force (Vtensor) is also shown. As 
was explained in the previous paper [35j, the properties of 
each minimum are characterized with their dominant 2p- 
2h configurations. At the lowest minimum point (a), the 
dominant 2p-2h configuration is given by (0s) 1 Q 2 (0p 1 / 2 )ioj 
where the subscripts 01 or 10 represent spin and isospin 
for the two-nucleon pair, respectively. This configura- 
tion provides the largest tensor force contribution with 
a spatial shrinkage of the Opi/% orbit and produces the 
high momentum component, similar to the results in 
Refs. (3(| UK]. At this energy minimum, (Vt cnsor ) is 
about three times as large as in the ordinary shell model 
case with bo Pl/2 /bo s —l. At the second minimum point 

(b), (0p3/ 2 )oi 2 (0pi/2)ou the Op shell pairing correlation is 
dominant. The results indicate that (a) and (b) represent 
the tensor correlation and the pairing one, respectively. 

As the second step, we superpose the above two wave 
functions at (a) and (b) in order to obtain the 9 Li wave 
function including the tensor and pairing correlations, si- 
multaneously. In the superposed 9 Li wave function, the 
dominant configurations are as follows with their prob- 
abilities Op-0/i with 82.9 %, (0p 3/2 )o 1 2 (0pi /2 )oi with 9.0 
% and (0s) ro 2 (0pi/2)io an d 7.2 %, respectively. From 
the two kinds of the preferred 2p-2h configurations, the 
obtained 9 Li wave function is found to possess both the 
tensor and pairing correlations. We show the contribu- 
tions of the Hamiltonian component in Eq. (TTJ) in Ta- 
ble U and the tensor force gives —20.7 MeV. The con- 
tribution from the pairing correlation is —3.0 MeV. We 
also show the results of the various r.m.s. radius R of 
matter(m), proton(p), neutron(n), charge(cft.) for 9 Li in 



5 




b 0p 1/2 I b 0s 

FIG. 1: Energy surface of 9 Li (solid line) with respect to the 
length parameters of Opi/2 orbit normalized by that of 0s 
orbit. The two minima indicated by (a) and (b) correspond 
to the states due to the tensor correlation and the paring 
correlation, respectively. The dashed line is the tensor force 
contribution. 



TABLE I: Contributions of the Hamiltonian to the Li ground 
state in units of MeV. 



Energy 


Kinetic 


Central Tensor 


LS 


Coulomb 


-45.3 


136.0 


-158.8 -20.7 


-3.6 


1.8 



Table [TTJ The recent observations of the charge radius is 
reproduced well in our wave function, where the radius 
of proton and neutron is considered. We present other 
properties for 9 Li in comparison with other Li isotopes 
later. 



B. Pauli-blocking effects in Li and Li 

Considering the properties of the configuration mixing 
of 9 Li, we discuss the Pauli-blocking effects in 10 Li and 
11 Li and their difference as shown in Fig. [2] For (a) in 
Fig. [2j the 9 Li ground state (GS) is strongly mixed, in 
addition to the Op-Oh state, with the 2p-2h states caused 
by the tensor and pairing correlations. 

Let us add one neutron to 9 Li for 10 Li. For (b) in 
Fig. [2j when a last neutron occupies the 0pi/ 2 -orbit for 
the p-state of 10 Li, the 2p-2h excitation of the pairing cor- 
relation in 9 Li are Pauli-blocked. The tensor correlation 
is also blocked partially, but not fully by Pauli-principle 
because the 0pi/2 orbit is not fully occupied by a last 
neutron. Accordingly, the correlation energy of 9 Li is 
partially lost inside 10 Li. Contrastingly, for (c) in Fig. [2j 
the Is state of 10 Li, the Pauli-blocking does not occur 
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TABLE II: Various r.m.s. radius of 9 Li given in unit of fm in 
TOSM and those of experiment. 





TOSM 


Experiment 


Rm 


2.31 


2.32±0.02f51] 






2.44±0.06[4] 


Rp 


2.10 


2.18±0.02[51] 


Rn 


2.41 


2.39±0.02f51] 


Rch 


2.23 


2.217±0.035[9] 






2.185±0.033[10] 



and 9 Li gains its correlation energy fully by the config- 
uration mixing with the 2p-2h excitations. Hence, the 
energy difference between p and s states of 10 Li becomes 
small to explain the inversion phenomenon [T^. |40| . 

For 11 Li, let us add two neutrons to 9 Li. The simi- 
lar blocking effect is expected for 11 Li, whose important 
properties were given in the previous paper (35j. For (d) 
in Fig. O when two neutrons occupy the 0pi/2-orbit, the 
2p-2h excitations of the tensor and pairing correlations in 
9 Li are Pauli-blocked. In particular, the blocking of the 
tensor correlation in n Li is expected to work stronger 
than the 10 Li case, due to the presence of the last two 
neutrons in the P\ji orbit. Accordingly, the correlation 
energy of 9 Li is lost inside 11 Li stronger than the 10 Li 
case. For (e) in Fig. [2J the (Is) 2 of two neutrons, the 
Pauli-blocking does not occur, similar to the Is state of 
10 Li. Hence, the relative energy distance between (Op) 2 
and (Is) 2 configurations of 11 Li becomes small to break 
the magicity in n Li. It is found that there is a differ- 
ence in the blocking effects between n Li and 10 Li. It is 
interesting to examine how this difference affects the s-p 
shell gap problem in these nuclei. For the pairing cor- 
relation, we already pointed out the different blocking 
effects between 10 Li and 11 Li in the previous study [H. 
We further consider the blocking effect in the dipole ex- 
cited states in n Li later, which is also different from the 
ground state case. In the previous paper (3^1, we exam- 
ined that the configuration mixing of the sd-shell for 9 Li 
would give a small influence on the blocking effect on the 
(Is) 2 configuration of n Li. 



'Li 



We solve 10 Li in a coupled 9 Li+n model, in which the 
9 Li-n folding potential is determined to reproduce the 
two-neutron separation energy S^n of n Li as 0.31 MeV. 
On this condition, we investigate the spectroscopic prop- 
erties of 10 Li. In Table Hill it is shown that using TOSM 
for 9 Li, the dual p-state resonances are obtained near 
the 9 Li+n threshold energy. The dual states come from 
the coupling of [(0p 3 / 2 ) w (0pi/ 2 )i/]i+/2+ , while the exper- 
imental uncertainty is still remaining including the spin 
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FIG. 2: Schematic illustration for the Pauli-blocking in Li 
and xl Li. Details are described in the text. 



assignment [26]. The 1 + state is predicted at a lower en- 
ergy than the 2 + state due to the attractive effect of the 
triplet-even 3 E channel in the pn interaction. 

For the s-wave states, their scattering lengths a s of 
the 9 Li+7T. system show negative values. In particular, 
the 2~ state shows a large negative value of a s , which 
is comparable to that of the nn system (—18.5 fm) with 
the 1 Sq channel, and indicates the existence of the virtual 
s-state near the 9 Li+n threshold energy. Therefore the 
inversion phenomenon in 10 Li is reasonably explained in 
the present model. The order of 2~ and 1~ also comes 
from the attractive 3 E component in the pn interaction. 

For comparison, we calculate 10 Li without the core ex- 
citations of 9 Li ("Inert Core"), namely, we adopt only 
the single Op-Oh configuration for 9 Li without the Pauli- 
blockings explained in Fig. ^d). In this case, we adjust 
the 8 parameter to be 0.066. In Table IIII1 the p-wave 
resonances are obtained just above the 9 Li+n threshold 
energy, and a s values show small positive values for both 
1~ and 2~ states, which means that the virtual s states 
are not located near the 9 Li+n threshold, and the s-p 
shell gap is large. These results mean that the Pauli- 
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TABLE III: The resonance energies E r and the decay widths 
r of the p-wave resonance states (1 + and 2 + states) in 10 Li 
measured from the 9 Li+n threshold are listed in unit of MeV. 
The scattering lengths a s of the s-states (1~ and 2~ states) 
are shown in unit of fm. We show here the two kinds of the 
results using TOSM and Inert Core for 9 Li. 





TOSM 


Inert Core 


(£ r ,r)(i+ 


) [MeV] 


(0.22, 0.09) 


(0.03, 0.005) 




) [MeV] 


(0.64, 0.45) 


(0.33, 0.20) 


a.(l") 


[fm] 


-5.6 


1.4 


a s (2") 


[fm] 


-17.4 


0.8 



blocking nicely describes the spectroscopic properties of 
10 Li. 

Here, we visualize the Pauli-blocking effect occurred in 
10 Li using the effective 9 Li-n potential. In Fig. [31 we de- 
pict the effective 9 Li-n potentials for the s-wave (V s ) and 
for the p-wave (V p ) of 10 Li, which are derived as follows; 
We use the inert description of the 9 Li core, namely, no 
Pauli-blocking from the tensor and pairing correlations 
and renormalize the blocking effect into V s and V p . For 
the s-wave, V s is given by the original folding potential 
V cn in Eq. JU) , because of no blocking from two correla- 
tions. On the other hand, for the p-wave, we reduce the 
strength of V cn to simulate the blocking effect by repro- 
ducing the resonance energy of the 1 + state (0.22 MeV) 
of 10 Li as obtained in TOSM shown in Table (ml This 
procedure is similarly done as explained in Refs. |19l.[40l]. 
In this manner, we obtain the state-dependent interac- 
tions with a deeper one for s-wave. The difference, Ap B , 
between V p and V s reflects the strength of the Pauli- 
blocking (PB): 

V p (r) ~ V cn (r) + A n PB (r), (15) 
V s (r) ~ V cn (r), (16) 
ApbM ~ V p (r)-V s (r). (17) 

Ap B depends on the relative distance r between the 9 Li 
core and a neutron, where superscript n indicates a last 
neutron in 10 Li. From Fig. [31 it is found that the Pauli- 
blocking acts repulsively for the p-wave, which explains 
the inversion phenomenon in 10 Li. 

We briefly mention the experimental situation of 10 Li. 
The experimental identification of the s-wave and p-wave 
states in the 9 Li+n threshold energy region are diffi- 
cult and the conclusive understanding of the 10 Li spec- 
troscopy is still an open problem [26| . In the recent obser- 
vation [26|. E r ~ 0.38 MeV with the decay width T ~ 0.2 
MeV is suggested as the p-wave resonance, which is rather 
close to our results for the 1 + state. For the s-wave state, 
a s = — 13 ~ — 24 fm is suggested, which agrees with our 
results of the 2~ state. It is highly desired that further 
experimental data are obtained for the low-lying states 
of 10 Li. 




r core-n [ fm ] 

FIG. 3: Effective 9 Li-n potentials, V s and V p obtained by 
renormalizing the core excitation in 10 Li, where the 9 Li-n po- 
tentials are transformed into equivalent ones using the inert 
9 Li core. V s is taken to be the same as the original folding 
potential V cn . 



D. n Li 

We solve 11 Li in a coupled 9 Li+n+n model and show 
the detailed properties of 11 Li in Table IIV1 the partial 
wave components P((nlj) 2 ) for halo neutrons, the vari- 
ous r.m.s. radius, the relative distance between halo neu- 
trons (Rn-n) an d the core-2n distance (R c -2n), and the 
expectation value of the opening angle between two neu- 
trons measured from the 9 Li core, respectively. The 
case using TOSM for 9 Li gives a large P((ls) 2 ) value, 
comparable to P ((pi/2) 2 ) for halo neutrons and a large 
matter radius R m for 11 Li, which are enough to explain 
the observations. The case of "Inert Core" gives small 
P((ls) 2 ) and small R m values, which disagree with the 
experiments. From the difference between two models, 
it is found that the tensor and pairing correlations in 
9 Li play important roles to break the magicity and make 
the halo structure of 11 Li, in addition to the s-p inver- 
sion phenomenon in 10 Li. As was shown in the previous 
study [35l]. the Pauli blocking effect from the tensor cor- 
relation is stronger than the pairing one to break the 
magicity of 11 Li. 

Here, similar to the 10 Li system, we show the Pauli- 
blocking effect in 11 Li using the effective 9 Li-n potential 
and discuss the difference of the blockings between 10 Li 
and n Li. In Fig. [H we depict the effective 9 Li-n po- 
tentials V s for the (Is) 2 configuration of n Li and V p 
for the (Op) 2 configuration, which are derived as follows; 
We use the inert description of the 9 Li core, namely, 

no Pauli-blocking from the tensor and pairing correla- 

2 

tions and renormalize this blocking effect into V s and 

2 2 
V p . For s- waves, V s is given by the original folding 

potential V cn in Eq. ([4]), which is also the same as V s 
in Fig. O For p- waves, we reduce the strength of V cn 
to simulate the blocking effects on the (Op) 2 configura- 
tion as follows; First, we take only the (Op) 2 configura- 
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TABLE IV: Ground state properties of 11 Li with S 2n = 0.31 
MeV, where two kinds of the 9 Li descriptions of TOSM and 
Inert Core are shown, respectively. Details are described in 
the text. 





TOSM 


Inert Core 


Expt. 


■P((Pl/2) 2 )) [%] 


42.7 


90.6 




(lsi/ 2 ) 2 


46.9 


4.3 


45±10[2l 


(p 3 / 2 ) 2 


2.5 


0.8 




(dsn) 2 


1.9 


1.3 




\ , - v o/ £ j 


4.1 


2.1 




ih/2) 2 


0.5 


0.2 




(A/2) 2 


0.6 


0.3 




Rm [fm] 


3.41 


2.99 


3.12±0.16[51] 








3.53±0.06f3] 








3.71±0.20f4] 


R p 


2.34 


2.24 


2.88±0.11[51J 


Rn 


3.73 


3.23 


3.21±0.17[51] 


Rch 


2.44 


2.34 


2.467±0.037[9] 








2.423±0.034[10] 


Rn-n 


7.33 


6.43 




Rc-2n 


5.69 


4.26 




[deg.] 


65.3 


73.1 





tion for x( nn ) in Eq. (|10[) using TOSM for 9 Li, which 
gives (E r , T) = (1.50, 0.65) MeV of the n Li ground state 
energy. Second, we adopt the inert 9 Li core and adjust 
V cn to reproduce the above resonance energy, which gives 
V p . The difference, A PB: between V p and V s reflects 
the strength of the Pauli-blocking for each last neutron 
in n Li: 

VP 2 (r) ~ V cn (r) + Ap i3 (r), (18) 
V s2 (r) ~ V cn (r), (19) 
Af B (r) ~ V*\r)-V'\r)> (20) 

where the superscript n 2 indicates the last two neutrons 
in n Li. In Fig.0] it is found that the Pauli-blocking acts 
repulsively on the (Op) 2 configuration of 11 Li and makes 
the (Op) 2 and (Is) 2 configurations placed at the degen- 
erated energies, so that they can couple by the pairing 
interaction. The energy difference between two config- 
urations is —0.1 MeV in 11 Li obtained in the previous 
study (3o|. In comparison with the 10 Li system shown in 

Fig. [31 it is found that the blocking effect A PB in 11 Li 
is stronger than A PB for 10 Li. This difference comes 
from the blocking of the tensor correlation in n Li being 
stronger than the 10 Li case, due to the presence of the 
last two neutrons in the Px/2 orbit, as was explained in 
Sec. ILUBl 

In Table [Vj we see the structure change of 9 Li inside 
11 Li from the relative energy difference AE between the 
9 Li core inside n Li and the isolated 9 Li ground state. 




Wn [ fm ] 

FIG. 4: Effective 9 Li-n potentials, V s2 and V p2 obtained by 
renormalizing the core excitation in 11 Li, where the 9 Li-n po- 
tentials are transformed into equivalent ones using the inert 
Li core. V s is taken to be the same as the original folding 
potential V cn . 



TABLE V: Internal energy difference between 9 Li inside the 
11 Li ground state and the isolated 9 Li ground state in units 
of MeV. 



AE 


Kinetic 


Central 


Tensor 


LS 


Coulomb 


2.0 


-0.8 


0.5 


2.7 


-0.3 


0.0004 



The differences of each component in the 9 Li Hamilto- 
nian are also shown. The positive value indicates the 
energy loss. Among the Hamiltonian of 9 Li, it is found 
that the energy loss of the tensor force is the largest. This 
is caused by the Pauli-blocking effect on the tensor cor- 
relation of 9 Li in n Li. Corresponding to this, the kinetic 
energy is gained due to the reduction of the high momen- 
tum component produced by the tensor force. Totally, 
the 9 Li cluster is excited in energy in the n Li ground 
state. 

The spatial correlations of the halo neutrons in 11 Li is 
interesting [l]| [H, [H, HiJ and shown in Fig. [SI We depict 
the density distribution of halo neutrons p nn (r, 0) in 11 Li 
as functions of 9 Li-n distance, r, and the opening angle 
between two neutrons, 0, with the following definition; 

/>oo 

Pnn{r,0) = / dr'p nn (r,r',8) (21) 
Jo 

Pnn(r,r',e) = 8ttW 2 smO 

x (y J ( n Li,r,r',6)\8(r - r")5(r' - r'") 
x 5(6-6')\y J ( n U,r",r'",9')}, (22) 

where the total wave function of n Li has three variables 
of the 9 Li-n distances with r and r' for each neutron 
and their opening angle 6 and in Eq. (f22|) . we integrate 
only the variables of the ket part. In the present TOSM 
case, (a), it is confirmed that the di- neutron type con- 
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'core 



(Color online) Two neutron correlation density 
n ,8) for n Li. (a) is the present calculation using 

respec- 



FIG. 5: 

Pnn (^*core — n 

TOSM of 9 Li and (b) is the Inert Core case of 9 Li 
tively. 



figuration ( a large r and a small 9) gives a maximum 
value of the density, although the density of neutrons is 
widely distributed. Contrastingly, the Inert Core case, 
(b) with a small s 2 component does not show much en- 
hancement of the di-neutron type configuration and the 
cigar type configuration (a small r and a large 9) coexists 
with it, which is similar to the other halo nuclei, such as 
6 He[52|. These two results indicate the role of the s 2 com- 
ponent on the formation of the di-neutron type configu- 
ration as follows: The s 2 component in the 11 Li ground 
state increases the amplitude of the tail region of two 
neutrons far from 9 Li, and these neutrons tend to come 
close to each other to gain the interaction energy between 
them. As a result, the di-neutron type configuration is 
enhanced, although the spatial distribution of neutrons 
are still wide. The spatial distribution of two neutrons 
also affects the opening angle 9, where the present TOSM 
case shows a smaller 9 value than the Inert Core one, as 
shown in Table [TV] 

The large s 2 component in n Li enhances the core-2n 
distance, which affects not only the matter radius but 
also the proton radius of 11 Li, because the proton radius 
of n Li consists of the core part and the core-2n distance 



> 
CD 



LU 
"O 

T3 




Present (S 2n =0.31 MeV) — 
(S 2n =0.38MeV) 



12 3 4 

Energy [MeV] 



FIG. 6: Calculated Coulomb breakup cross section measured 
from the 9 Li+n+n threshold energy in comparison to the 
experiments^. There are two lines with different S2n for 
the 11 Li ground state. 



{Rc-2n) in the following form; 



R^U) = i? 2 ( 9 Li) 
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R 



c-2n- 



(23) 



The latter term represents the recoil effect of the cm. 
motion from the halo neutrons. Experimentally, the 
charge radius of n Li was measured recently and its value 
is enhanced from the one of 9 Li[|[lo|, as shown in Table 
IV! The present wave functions provide reasonable value 
for 11 Li. This enhancement is mainly caused by the large 
distance between 9 Li and the paired 2n. The Inert Core 
case gives the small charge radii for 11 Li due to the small 
distance of R c -2n, since the s 2 component is very small. 

We show the three-body Coulomb breakup strength of 
11 Li from the ground state into the 9 Li+n+n system and 
compare it with the new data by the RIKEN group @ in 
Fig. [6l We use the Green's function method combined 
with the complex scaling method to calculate the three- 
body breakup strength[35l. 1551 ] using the dipole strength 
and the equivalent photon method, where experimental 
energy resolution is taken into account [8J. In Fig. [6l it 
is found that the present model reproduces the experi- 
ment, in particular, for the low energy peak observed at 
around 0.5 MeV. For reference, we show the other re- 
sults using a different two-neutron separation energy Sa n 
of the 11 Li ground state as 0.38 MeV suggested by the 
recent experiment (56j . Here we adjust the parameter 5 
in the 9 Li-n potential. In that case, the matter radius 
of n Li is 3.28 fm. In the strength, only the magnitude 
at the peak is slightly reduced from the present results 
using S , 2„=0.31 MeV due to the threshold effect from the 
halo structure in the ground state [13, H^|. Seeing more 
closely, our results seem to underestimate the cross sec- 
tion at E > 1 MeV, which is commonly seen in the results 
with different S^n- 

Several three-body calculations of 11 Li predicts the 
dipole resonances in the low excitation energy region 
with the decay widths, relatively smaller than the res- 
onance energies 21, 59]. In fact, when we assume the 
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TABLE VI: Q moments of 9 Li and n Li in units of e fm 2 
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FIG. 7: Schematic illustration for the Pauli-blocking in the 
dipole states in Li corresponding to the three kinds of the 
9 Li configurations. Details are described in the text. 
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FIG. 8: Effect of the final state interactions on the El tran- 
sition strength of Li. 



inert 9 Li core and take a phenomenological 9 Li-n poten- 
tial with a deep s-wave potential, the dipole resonances 
of 1/2+, 3/2+ and 5/2+ are obtained 35]. Contrastingly, 
in the present model of 11 Li with the tensor and pair- 
ing correlations, we cannot find any three-body dipole 
resonances 35] . Here, we discuss why the three-body 
dipole resonances are not obtained in the present model. 
To do this, we consider the Pauli-blocking effect on the 
dipole excited states in the analogy of the ground state as 
was explained in Sec. IIII Bl In Fig. [3 we show the dom- 
inant shell model type configuration for the dipole ex- 
cited states, in which the last two neutrons occupy lsj/2 
and 0pi/2 orbits, respectively. In this configuration, due 
to the Pauli-principlc, the blocking effect for the tensor 
part becomes weaker than that in the p 2 configuration 
of the ground state, shown in Fig. [DJd). This weakness 
of the blocking reduces the repulsive effect on the p-orbit 
and makes the energy gap between 0pi/ 2 and lsi/ 2 orbits 
wider to get close to an usual order of the shell model or- 
bits. In this situation, the (0pi/2)(lsi/2) configuration 
loses the energy, which means that the dipole resonances 
have relatively higher excitation energies than the phe- 
nomenological potential model. Furthermore, since the 
s-wave neutron is easy to decay, the dipole resonance has 
a large decay width when the resonance energy is rela- 
tively high, which makes one difficult to find out its loca- 
tion in energy spectra to confirm the physical signature, 
such as the peak in the strength function [551] . 





9 Li 


"Li 


TOSM 


-2.65 


-2.80 


AMD [60] 


-2.66 


-2.94 


SVM[18] 


-2.74 


-3.71 


Exp eriment fill 


-2.74 ±0.10 


3.12 ±0.45 (]<2|) 


Exp eriment [61] 


-3.06 ±0.02 





TABLE VII: n moments of 9 Li and 11 Li in units of fj,N. 





9 Li 


n Li 


TOSM 


3.69 


3.77 


AMD [60] 


3.42 


3.76 


SVM[18] 


3.43 


3.23 


Exp eriment [11] 


3.44 


3.67 



We further try to understand the physical origin of the 
low energy peak in the dipole strength around 0.5 MeV 
from the 9 Li±n±ra threshold energy. Here, we investigate 
the effect of the final state correlations induced by the in- 
teractions between the 9 Li±n±n system. Using TOSM 
for 9 Li, in Fig. [5J we show the other three strength func- 
tions in addition to the original full calculation, in which 
either the core-n or nn interactions is cut off, respectively, 
or both of them are cut off, namely, the plane wave treat- 
ment of the dipole excited states. In these calculations, 
we do not change the ground state properties and do not 
take care of the experimental energy resolution of the 
strength. The results show that the plane wave gives the 
small peak at 0.5 MeV in comparison with the present 
model, which indicates the importance of the final state 
correlations. Among the final state correlations, the core- 
n and nn interactions give comparable contributions and 
both effects are necessary to reproduce the low energy 
peak in the TOSM case. This fact means that both cor- 
relations of core-n and nn are important to explain the 
breakup mechanism of 11 Li into the 9 Li±n±n system. 
The effects of the p-wave resonances and the virtual s- 
states of 10 Li, and of the di-neutron type configuration in 
the final states are interesting and would be investigated 
elsewhere. 



E. Electromagnetic properties of Li isotopes 

We show the Q and p moments of the 9 Li and 11 Li 
in Tables IVII and IVII1 where only the absolute value of 
the Q moment of n Li is reported in the experiments [TTI] . 
The present model describes reasonably those values for 
9 Li and 11 Li. For Q moments, the values of 9 Li and n Li 
do not differ so much to each other. This result is similar 
to that of AMD[60[ and different from SVM based on 
the multi-cluster model[l^|. Here, similar to the charge 
radius, we discuss the recoil effect in the Q moment of 
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11 Li by expanding its operator Q( 11 Li) into the core part 
(Q( 9 Li)), the last neutron part and their coupling part 
as 



Q( n Li) = Q(»U) + ] /^3e(J^j y 2 o(Rc-2n) 

_ 87ry|[Oi( 9 Li),y 1 (R c . 2 „)] 20! (24) 
Oi m ( 9 Li) = e y^i), (25) 

ifEproton 

where 3^ m (r) = r Y; m (f) and {a^} are the internal coor- 
dinates of protons in 9 Li. In our wave function of 11 Li, 
last two neutrons almost form the + state with the prob- 
ability of around 99%. In that case, the Q moment for the 
relative motion part of 9 Li-2n(c-2n) having the relative 
coordinate R c - 2 n almost vanishes because of the non-zero 
rank properties of the Q moment operator. This means 
that the recoil effect from the clusterization is negligible. 
Therefore, the Q moment of n Li is caused mainly by the 
9 Li core part inside n Li. In our wave function, the spa- 
tial properties of the proton part of 9 Li inside n Li do not 
change so much. Hence, the Q moment of Li is similar 
to the value of the isolated 9 Li. Small enhancement from 
9 Li to n Li mainly comes from the lacking of the high 
momentum component of the tensor correlation due to 
the Pauli-blocking in 11 Li, which extends the radial wave 
function of n Li. The experimental information of the 
Q moment is important to understand the structure of 
11 Li. It is highly desired that further experimental data 
are available for the Q moment of n Li. 

For the fj, moment, the value in n Li is almost the 
Schmidt value of 3.79 of the 0p 3 / 2 proton. In 9 Li, 
the P1/2 proton is slightly mixed, which decrease the \x 
moment. In 9 Li, this P1/2 proton is excited from the 
0p 3 / 2 orbit in a pair with the P1/2 neutron and the 2p- 
2h excitation is Pauli-blocked in 11 Li due to the addi- 
tional neutrons. As a result, the excitation of P1/2 is 
suppressed, which makes the fi moment of 11 Li close to 
the Schmidt value of the p 3 / 2 orbit. The tendency of 
the increase of the fi moment from 9 Li to 11 Li can be 
obtained also in the shell model analysis using various 
effective interactions [62j. 



IV. SUMMARY 

In summary, we have investigated the structures of 
9,10, Hjj systematically based on the extended three-body 
model, in which the tensor and pairing correlations are 
explicitly considered for the 9 Li part. It is found that the 
two correlations in 9 Li play an essential role to explain 
the anomalous structures of 10 ' n Li. The obtained results 
clearly demonstrate that the inert core assumption of 9 Li 
is not reliable. For 10,11 Li, the Pauli-blocking induced 
by the tensor and pairing correlations in 9 Li naturally 
explains not only the inversion phenomena of 10 Li, but 
also the breaking of magicity and the halo formation in 
n Li consistently. It is pointed out that the strengths of 
the Pauli-blocking between 10 Li and 11 Li including their 
excited states are different. The interplay between the 
tensor correlation of 9 Li and the halo neutrons mainly 
gives rise to this difference and due to this difference, the 
blocking acts strongly for the n Li ground state. We have 
further investigated the various physical quantities of Li 
isotopes. In n Li, it is shown that the development of the 
halo structure induces the di-neutron type configuration, 
and the large matter and charge radius. The electro- 
magnetic properties of 9 Li and n Li are reproduced in 
the present model. The present model explains further 
the recent observation of the Coulomb breakup strength 
of n Li. The low-energy peak in the Coulomb breakup 
strength originates from the final state interactions. 
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